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Estimation of halo ellipticity using spin-3 flexion 
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ABSTRACT 

Estimating the ellipticity of dark matter haloes at the galaxy or galaxy-cluster scale can 
provide important constraints on the formation of galaxies or clusters, as well as on the 
nature of dark matter. We show in this paper that the spin-3 gravitational flexion can add 
useful information on the ellipticity of lensing haloes. We introduce a general formalism 
to decompose fields with arbitrary spin into radial and tangential components. The ratio 
of the tangential and radial flexion components directly estimates the lens ellipticity. We 
point out that any centroid offset will significantly bias our estimate, which on the other 
hand can be used to determine the centre of the lens halo. 
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1 INTRODUCTION 

The structure of dark-matter haloes hosting clusters and galax- 
ies is of substantial interests for cosmology. Halo properties con- 
sistently predicted by A'-body simulations of structure forma- 
tion in Cold Dark Matter (CDM) include a highly non-spherical 
structure, which is wel l fitted by a triaxial body with a charac - 
teristic density profile iJing & Sutdl2002l ; lspringel et al.ll2004 . 
This property is generic for the nature of the dark matter 
as well as for the forma tion of galaxies and clusters (B ullock 
l2002r,lAllgood et al.l2006l ; lBett et al.l201ol ; IVera-Ciro et al.l201 1| ; 
Wang et afll20ld . l201ll) . For example, major mergers of dark- 
matter haloes may play an important rol e in forming their shapes 
Ivan Haarlem & van de Wevaaertjl993l) . Testing halo ellipticities 
for their agreement with the cold-dark matter prediction may thus 
provide a potentially powerful constraint on the standard sce- 
nario for structure formation from galaxies to galaxy clusters. 

Arguably the most promising tool to study the mass dis- 
tribution of galaxy and cluster haloes is gravitational lensing 
iBartelmann & Schneiderll200ll) , as it probes the mass distribu- 
tion independent of the nature of the matter (luminous or dark) or 
its dynamical state. Weak gravitational lensing is routinely bein- 
used for cluster m ass reconstructions (e.g. Iciowe et alJlioQi 
Brad ac et ail I2008I) and for measuring the ellipticity of dark- 
matter haloes j Mandelba um etal.l l2006; Corless & Kind l2007l ; 
iDeb et al.l2O10l ). For example, a mean ellipticity of 0.46 has been 
obtained f or a sample of 25 galaxy-cluster haloes based on Sub- 
aru data jOguri et all 120101 ), roughly consistent with the CDM 
prediction. 

Hitherto, the shear has mainly been used for weak-lensing 
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studies, and the gravitational magnification in a few cases. Both 
are determined by the projected gravitational tidal field, i.e. by 
linear combinations of second-order derivatives of the effective 
lensing potential. Measurements of third-order derivatives are 
now gradually coming within reach. They can conveniently be 
combined into quantities called the first and the second flexion, 
or T and Q flexion, respectively. More on their properties will 
follow below. 

The two flexions can be introduced as derivatives of 
either the surface mass density or the gravitational shear. 
They respond to smaller-scale variations in the projected 
mass distri bution than the s hear a nd the surface-mass den- 
sity jGoldberg & Nataraianl l2002r , iGoldberg & Baconl 120051 : 
iBacon et al.ll2006l) . S everal techniques have been proposed to 
measure flexion (e.g. lOkura et al.ll2007l ; IVeTander et al.ll201 1| ; 
ICain et al]|201ll) , but an unambiguous measurement has yet to 
be demonstrated. It has recently been shown that a reliable es- 
timate of the flexions requires an accurate model for the shear, 
which must be taken in to account in any unbiased flexion esti- 
mator jViolaetal.1120121) . 

The intrinsic noise of measured flexion com ponents is 
still poorly understood (Goldberg & Leonard |2007|) . However, 
it has been suggested that flexi on can signif icantly improve 
the estimation of hal o ellipticity lEr et al.ll2012l) . In particular, 
lEr & Schneider! |201 1|) have introduced the ratio of the tangential 
and radial components of the T flexion to estimate the halo el- 
lipticity. Since this approach is independent of the lens strength, 
and thus also independent of the mass of the lens, its redshift 
and even the halo-centric distance, one should be able to use 
this ratio of T flexion c omponents to estimate the halo ellipticity 
as a function of radius lEr et alJlioi 

Here, we extend this earlier study in three ways, motivated 
by the following considerations: First, the T flexion is a quantity 
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with the transformation properties of a vector. It is thus straight- 
forward to define its radial and tangential components. The Q 
flexion has a three-fold rotational symmetry, which renders its 
projection on tangential or radial directions less transparent. 
We propose here a general formalism for decomposing lensing 
quantities of arbitrary order and spin weight into components of 
different orientations, which may also help clarifying the proper- 
ties of the T and Q flexions. Second, the T flexion leads to a 
centroid shift in images which is probably hard to measure re- 
liably. The Q flexion with its unique, triangular distortion pattern 
should be considerably easier to access. It is therefore important 
to exploit in particular the signal of the Q flexion, if flexion can be 
measured at all. Third, since the signal is expected to be weak, 
the combination of all accessible signals is mandatory, calling for 
the Q flexion to be included in measurements of halo shapes. 

In Sect. 2, we introduce our formalism for decomposing 
spin-weighted fields of arbitrary spin into tangential and radial 
components, and apply it specifically to the gravitational flexion. 
We show the results of numerical tests in Sect. 3. We discuss 
these results and summarise our conclusions in Sect. 4. Where 
necessary throughout this paper, we adopt a ACDM standard 
cosmological model with n A = 0.75, Q m = 0.25, and a Hubble 
constant of H = 73 km s~' Mpc 1 . 



imaginary parts by (T,ff)i and (T,Q) 2 , respectively. In terms of 
the lensing potential, we have 

T = Ti + iT 2 = - (<9fy + d&ijr) + l - (djd^ + d\ij,) (5) 



and 



g = g x + \g 2 = x - (dfy _ 33i3^) + i (3d 2 l a 2 i/> - a^) . (6) 



2.2 Spin decomposition of tensors on the sphere 

We proceed by introducing a general formalism for decompos- 
ing arbitrary lensing distortions, or more generally fields on 
the sphere with arbitrary spin, into fields of defined rotational 
symmetry and projecting them onto arbitrary dire ctions on the 
sky. F or reference, we recommend Appendix A of ICastro et alj 
12005) and the literature cited therein. Let e 12 be the conven- 
tional coordinate basis vectors on the sphere, 



e 2 = sin 88, 



(7) 



and 8'~ their dual basis vectors, 



g 1 =d9 , 6 2 = sinSd^ . (8) 
We introduce the conventional helicity basis e± and its dual basis 



2 FORMALISM 
2.1 Preliminaries 

The fundamentals of g r avitati o nal lensing can b e found in 
iBartelmann & Schneider! l200ll) ; lBacon et al.l l2006ft . For its el- 
egance and brevity, we shall use the complex notation for shear 
and flexion. The thin-lens approximation is adopted, implying 
that the lensing mass distribution can be projected onto the lens 
plane perpendicular to the line-of-sight. We introduce angular 
coordinates 6 with respect to the line-of-sight. The lensing con- 
vergence, i.e., the dimensionless projected surface-mass den- 
sity, can be written as k0) = E(f?)/E cr , where 2,0) is the projected 
surface-mass density and S ct is its critical value 
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depending on the angular-diameter distances £> s , Ai and D ds 
from the observer to the source, the observer to the lens, and 
the lens to the source, respectively. All lensing quantities can be 
derived from the effective lensing potential ifi, 



- f d 2 eV(0')ln|0 - ff\ 

7T Jr2 



(2) 



To lowest order, image distortions caused by gravitational lens- 
ing are described by the complex shear 



y = 7i + m = -z [dli// - djif/) + ididzi// 



(3) 



The shear transforms a hypothetical round source into an el- 
liptical image. The T and Q flexions can be introduced as the 
complex derivatives 



T = (di+id z )K, 
Q = (<9i+id 2 )y. 



(4) 



The flexions are thus combinations of third-order derivatives of 
the effective lensing potential iff. We shall denote their real and 



e ± = —= Oi ± ie 2 ) , e ± = — (e l + ie 2 ) 

V2 V2 v ; 



(9) 



The basis vectors e+ have spin ±1. From the definition of these 
bases, the relations 

6 1 (e ± )=-Lff k (e l ), e 2 (e ± )=^(e 2 ) (10) 
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follow easily since ff(ej) = {?. for i, j = 1,2 (or (', j = +, -). Let now 
T be an arbitrary rank-r tensor field, then this tensor applied to 
r copies of e±, 



± ,-t = T(e ± ,...,e ± ) 



(11) 



is a function with spin ±r. Rank-r tensors with defined spin prop- 
erties can thus be constructed as 



f = (J) 6 + ®...®6 + + (_,.?) er 



(12) 



Fully symmetric tensors, which we are primarily concerned 
with, can be expanded into the ordered sum 



'■®...®9 2 , (13) 



where an exponent k on an index denotes k repetitions of this 
index. The spin-±r functions ±r t are then 



^=2^g(3 r "- i2i(±i)t - 



(14) 



Of particular interest for us are second- and third-rank, symmet- 
ric tensors. For r = 2, we find 



±2?= ^(T n -T 22 )±\Tn 



(15) 



while r = 3 gives 



±3 f = j!?! [(Tin - 3r 122 )± i(3r 112 -r 222 )] . (16) 

Quite obviously, Eqs. (15) and (TB) reproduce the shear and the 
^-flexion, respectively, and their complex conjugates. 
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At any arbitrary point on the sphere, we now introduce a 
polar angle <p with respect to the local e\ direction and define a 
tangential and a radial unit vector, e T and e R , respectively, by 

e R = R(<p)e, = —R(ip)(e+ + e_) , 
V2 

e T = R(<p)e2 = —^=R(<p)(e + - e.) , (17) 

where R(<p) is a matrix representing a (two-dimensional) rotation 
by an angle <p. Since e ± are defined to be eigenvectors of R(jp) 
with eigenvalues e TV , 

«* = ^ (e-'V + e'^_) , e r = ( e -^ e+ - e '^_) , (18) 

the dual helicity basis vectors 9 ± applied to these vectors give 



V2 V2 



(19) 



The radial and tangential components of the spin-±r tensor field 
f can now be defined by the application of f to r copies of the 
radial and tangential unit vectors e R and e T , respectively, 

T R = T(e R ,...,e R ), T T = T(e T ,...,e T ). (20) 

With the expansion Eq. (T2), this gives 

f R = ^ + M)e^] , 

fr = ^ [i r G?)e^ + (-i) r (-,-0e^] . (21) 

If we insert the expressions for ±r t for symmetric tensors from 
Eq. fi"4l here, we finally obtain 

f - = ^Z(l) r ^ i,!l ( iV, 1 < 22 > 



for the radial component and 



(23) 



for the tangential component. Quite obviously, for even rank r, 
these components are related by f T = i'f R . For a symmetric, 
second-rank tensor, these general Eqs. (22) and (23} specialise 
to 

t R = - (Tn - T 22 ) cos 2 V + T 12 sin 2tp , f T = -f R . (24) 



2.3 Application to gravitational lensing 

We now proceed to apply this general formalism to second- and 
third-order partial derivatives of the effective lensing potential if/. 
With the second-order partial derivatives of the lensing potential 
if), we can define the rank-2 tensor 



(25) 



from which we can project out functions with spin and spins ±2 
by 

0l // 2 > = <f, {2 \e + ,e_), ±2 ^=^ 2 '( e± , e± ). (26) 

Straightforward calculation shows that the spin-0 function is the 
lensing convergence, 



0<A (2) =*. 



(27) 



while the spin-±2 functions are related to the shear by 

2 l/r <2) = y = y l +iy 2 , _ 2 | A <2> = 7* = 7l ~ iy 2 • 



(28) 



Accordingly, we define a rank-2 shear tensor y from the spin-±2 
components of the tensor if/ (2) introduced in (T2) above, 



y = y0 + ®Q + +y'9 ®8 



(29) 



The usual tangential and cross components of the shear, 
y T and y x , are then straightforwardly introduced as the results 
of the tensor y applied to (e T , e T ) and (e T , e R ), 



y T 



: %e T , e T ) = -% ( r e- 2i ^) , 7x = y{e T , e R ) = 3 (ye^) 



(30) 



which reproduce the common definitions 

y T = -y\ cos 2<f — y 2 sin 2ip , y x = -y 2 cos 2<p + y\ sin 2ip . (31 ) 

Applying the general formalism to the flexion now, we begin 
with the rank-3 tensor 



f// 3) = didjd k if/8' <S> 6^ ® 6»* , 
and extract its spin-1 component by the operation 

Il /r (3 ' = l /r (3) (e + ,e + ,e-)- 



(32) 



(33) 



A brief calculation shows that the result is related to the In- 
flexion by 



^ (3) = 



V2 



(34) 



Since the y-flexion is a vector, its radial and tangential compo- 
nents are trivially given by 

T R = f\ cos <p + Ti sin ip , Tt = fx sin (p - !F 2 cos </> . (35) 

Specialising the general Eqs. (22) and (23) to r = 3, and multi- 
plying by a factor of two for convenience, we find the radial and 
the tangential components of the ^-flexion, 

Sr = @i cos 3tp + §2 sin3y> , 

@t = @\ sin3i/> - @2 cos 3ip . (36) 

2.4 Singular isothermal elliptical halo 

Having laid out a general formalism for defining the tangential 
and radial components of a symmetric rank-r tensor field on the 
sphere, we now proceed to calculate the radial and tangential 
components of the ^-flexion for a dark-matter halo modelled as 
a Singular Isothermal Ellipsoid (SIE) with ellipticity e. The halo 
ellipticity is defined by e = (d a - 9 b )/(d a + 9 b ), where 8 a and b 
axe the major and minor axes, respectively. The dimensionless 
surface mass density for a SIE halo is 



0e(1 - e) 
2p 



(37) 



where p = yjpd 2 + 2 is a radial coordinate which is constant 
on elliptical contours, and / = 9 b /9 a is the axis ratio. Then from 
our result (36) for the two ^-flexion components, we have 

3 6>e(1 -e) 
2 



Sr 



9p 

-(l- e) (l-/^ 



(38) 
(39) 
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where 8 = ^6] + 9\ is the ordinary circular radius. The ratio be- 
tween the tangential and the radial components of the ^-flexion 
is similar to that of the T-flexion, 



Qt 2e sin 2ip Tt 

~Gr = 3(1 -26 cos 2 V + e 2 ) = Wr 



£r . (40) 



except that the ratio for the ^-flexion is one third times that for 
the T-flexion. In case of the SIE halo mass distribution, Q R is 
three times larger than T R . However, the tangential components 
of both the T- and ^-flexions are only affected by the ellipticity of 
lensing mass distribution, i.e. they are equal in magnitude. Thus 
the ratio between the components of the ^-flexion is one third 
times that of the T-flexion. Thus, if we define the ratio for the 
^-flexion as q = 3(0t/@r), an identical approach as for the in- 
flexion can be employed also for the ^-flexion. In particular, the 
estim ators {g) and ( g 1 ) for the halo ellipticity are then the same 
as in lEretal'llioTll) . 

Ignoring noise, the ratio q for the T- and ^-flexions will be 
identical. Then, the ratio for the ^-flexion would be redundant 
with the ratio from the T-flexion and would thus not provide ex- 
tra information. With noise, however, we expect two sets of data 
from the T- and ^-flexions with different noise properties. In par- 
ticular, since the ^-flexion causes a triangular-shaped distortion 
of images, its signal-to-noise ratio may be substantially higher 
than for the T-flexion which causes essentially a centroid shift 
of the background galaxy image. 

If the noise contributions to the T- and ^-flexion measure- 
ments can be considered independent, we can simply add the 
signal from T and Q, obtain twice the amount of data and 
achieve a lower sample variance. In reality, as pointed out by 
Viola et al.(2011), the shear needs to be known for an unbi- 
ased measurement of both types of flexion. Uncertainty in the 
shear may thus bias the flexion and may couple the measure- 
ment noise in the different types of flexion. The noise contribu- 
tions to the two types of flexion may thus be correlated, compli- 
cating the joint estimation of halo ellipticities from both types of 
flexion. However, lacking any detailed knowledge on the flexion 
noise, a simplified approach may be justified. 

In absence of noise, T T = Qt and 3T R = Q R . The data 
points will be particularly noisy if the tangential or the radial com- 
ponents of the T and § flexions are very different, i .e. if the noise 
dominates the measurement. Thus, we only use the simulated 
data if T T and Q T as well as W R and Q R are not very different, 
e.g. 0.5 < TtIQt < 2. This may be appropriate without detailed 
knowledge of the flexion noise. We do not aim at reducing the 
noise, but only disregard the most noisy data. The information in 
the most noisy data will thus not be used, reducing the amount 
of independent data. 



3 NUMERICAL TESTS 

In order to test the estimator for the halo ellipticity from the Q 
flexion, we use an SIE halo (Eq.[37) to generate mock data. We 
choose a lensing halo with an Einstein radius of 9 E = 5" and 
place it at redshift of z = 0.4. Background source galaxies are 
randomly distributed within a field sized V x 1' behind the lens. 
The redshifts of the source galaxies are drawn from a Gamma 
distribution with za = 1/3, 



which peaks at z = 2/3 and has a mean redshift of (z) = 3zo = 1. 

We create mock data sets with different values for the halo 
ellipticity, 0.1 < 6 < 0.6. For each ellipticity, we produce 50 re- 
alizations, with a spatial number density of noisy background 
galaxies of 40arcmin~ 2 . The average estimate over 50 realiza- 
tions is used as expected estimation, and the standard deviation 
is used as the errorbar. Lacking a reliable noise measurement, 
we adopt a simple model to include noise into our simulated 
data, 



0° 



Q\ + n s i + K&2 + n g2 ) . 



(42) 



The noise contributions n gl and n g2 are drawn from Gaussian 
distributions with identical variance, cr g i = a gl = 0.04arcsec~' 
for each component. Synthetic data for the T-flexion are also 
generated with Gaussian intrinsic noise characterized by th e 
variance 071 = 0.03 arcsec -1 = 072 1 Gold berg & Leonarj|2007l) . 

Flexion data are discarded if they are taken closer to the 
halo center than 4" or farther from the halo center than 14", 
where the flexion signal will be very small and the noise is 
expect ed to dominate . Sinc e large flexions cannot be mea- 
sured ISchneider & Erl 12008ft , data with \@\ > 1.5arcsec~' or 
|T| > 0.5arcsec~' are also disc arded. Moreover, due to the 
upper bound on the flexion ratio lEr & Schneider! 1201 ll) , ratios 
q > 3 are discarded as well. This amounts to assuming that the 
halo ellipticity does not exceed e = 0.7. We have approximately 
10 flexion data points for each realization on average. The flex- 
ion ratio q is calculated for each simulated flexion measurement. 
The mean flexion ratio is obtained by q = (1/AOZ/Ii 9u where 
N is the number of images for each realization. We use the esti- 
mator 



exp[nq/2] — 1 
exp[nq/2] + 1 



(43) 



P(z) 



—-exp(-z/z ) 

2Zn 



(41) 



for the halo ellipticity and ignore any redshift information in our 
tests. 

We perform two series of tests with our estimator and show 
our results in Fig. [TJ In one series, the noise of T and Q in the 
data are generated independently (left panel). In another series 
of data, we use correlated noise levels, namely n gl = (4/3)n fl 
and n g2 = (4/3)«/ 2 (right panel). In both panels, the solid line 
shows the input value. The stars show the result using the T- 
flexion, the open squares show the results using the ^-flexion, 
and the filled circles show the result using both the T- and Q- 
flexions. Note that very discrepant data points are discarded as 
described in the previous section in the combining analysis us- 
ing T and @. We can see that the estimate using the ^-flexion 
is significant larger than the estimate using the T-flexion in both 
panels. The reason for is that the ^-flexion ratio signal is one 
third of that of T-flexion. The noise thus affects the ^-flexion ra- 
tio more strongly. Another reason is that the noise we assume 
for ^-flexion is larger than the noise of T-flexion. In both panels, 
one can see that the estimations using two types of flexion pro- 
vide better results, i.e. lower bias, since very noisy data are dis- 
carded. On the other hand, the error bars of different estimations 
have no significant difference. In the combined estimation using 
T and Q, since we disregard discrepant measurements, the to- 
tal number of data is not twice as that of using only one flexion. 
Thus the error bars are not significantly smaller than those using 
only one type of flexion. 
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Figure 1. Comparison of the ellipticity estimator with the input ellipticity (solid line) for an SIE halo. The ellipticity is estimated from 50 realizations. The 
error bars represent the standard deviations. In the left (right) panels, the results obtained with uncorrelated (same) noise are shown. The red stars (blue squares, 
purple circulars) represent the results using the T^-fiexion (^-flexion, or both). 



3.1 Centroid offset noise 

Since the tangential and radial flexion components are calcu- 
lated with respect to the mass center of the lens, a centroid offset 
of the lens will affect the determination of the halo ellipticity. This 
bias has been studied and shown to be small using the F-flexion 
|Er & Schneiderll201 it) . Here, we test the effect of centroid off- 
sets on halo-ellipticity estimates from the ^-flexion. We simulate 
this effect by calculating the flexion ratio with respect to an as- 
sumed, perturbed centroid position d + 66, where 9 is the true 
center of the lens and 66 is the offset. We apply separate and 
independent offsets along the major and minor axes of the lens. 
To isolate this effect, we generate 200 realisations of halo fields 
without intrinsic noise on the source galaxies for a lens with el- 
lipticity 6 = 0.3. The results are shown in Fig. [2] The left (right) 
panel shows the estimate for the centroid shift applied along the 
major (minor) axis of the elliptical halo alone. The solid squares 
(crosses) present the estimate using the T- {§-) flexions. We 
can see that the bias due to the centroid offset is significant us- 
ing the ^-flexion. A similar result is obtained with noisy data. 
Moreover, the bias due to the centroid offset dominates the ef- 
fect of the intrinsic noise for the ^-flexion. Thus, considering the 
centroid offset, using the ^-flexion is more problematic than us- 
ing ^-flexion, but the ^-flexion may provide a better tracer for 
the center of lens halo. 

We can identify the center of the lens halo using the tangen- 
tial components of both types of flexion, T T and Q T , since the 
radial flexions are less affected by the centroid offset. We thus 
calculate the mean tangential flexion over all the background 
galaxies with respect to all the candidates of lens center. Here 
we choose 10000 points on a 6 x 6 arecsec 2 grid which is cen- 
tered at the true center of the lens. The center is then deter- 
mined at the position where we obtain the minimum mean tan- 
gential flexion. In Fig.|3]we show our estimated halo center off- 
set as a histogram using 200 noisy realizations for a halo with 
e = 0.3. Most of our estimates of the lens center have errors 
smaller than one arc second. We need to point out that this is 
valid for a regular SIE halo without any substructures. In case of 
non-relaxed cluster haloes, the mass center is often ill-defined 
and more difficult to identify. 
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Figure 3. Histogram of the estimated halo center offset for an SIE halo 
with e = 0.3 from 200 realizations. 

4 SUMMARY 

In this paper, we have introduced a general formalism for decom- 
posing lensing quantities of arbitrary order and spin into compo- 
nents of different orientation. We have applied the formalism to 
the spin-3 galaxy-galaxy flexion with elliptical mass distributions. 
We have derived the ratio of tangential to radial flexion and given 
the relation to the halo ellipticity. This ratio is similar to that for 
the spin-1 flexion, i.e. the ratio depends solely on the elliptic- 
ity and orientation of the lens mass distribution. Thus the spin-3 
flexion can provide extra information in the estimation. The noise 
behavior of the two types of flexion is still unclear, i.e. how large 
is the intrinsic noise and the noise correlation are for and across 
the T and Q flexions. We can only apply a simple data selection 
by comparing two types of radial (tangential) flexion and discard- 
ing the very noisy data. In the numerical tests, dropping the very 
noisy data can reduce the bias of our estimations. However, this 
method reduces the amount of data and loses information. Any 
specific noise cross-correlation between the two types of flexion 
can be employed in the future noise reduction, and a smaller 
bias will be expected with further noise analysis. 

We have found that unlike the spin-1 flexion, a misidentifi- 
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Figure 2. Bias due to a centroid offset of the lens. The solid line is the input value. The left (right) panel shows the result for a centroid shift applied to the 
major (minor) axis only. The solid squares (crosses) are the estimates using the T- (Q-) flexions. 



cation of the mass centroid can introduce a significant bias in 
our estimates using spin-3 flexion. We thus propose to use the 
mean tangential flexion to identify the mass center of the lensing 
halo. In the test using a regular SIE halo without subhaloes, the 
majority of the estimated centroid offsets using our method are 
smaller than one arc second. In the more realistic case, e.g. non- 
relaxed halo with substructures, further study is needed, possi- 
bly to be combined with other types of data, such as the X-ray 
surface brightness. 

The shear signal can also provide constraints on the halo 
properties and has a larger signal amplitude. However, shear 
should suffer more from the intrinsic shape noise of the back- 
ground galaxies. It can be estimated with the simple case of a 
singular isothermal halo profile with an Einstein radius t9 E . Let 
e E = 1 arc second, for example. The shear signal at = 59 E 
is 0.1, while the flexion signal is 0.02/arcsec. Since the intrin- 
sic standard deviation of the ellipticity is typically found to be 
cr e ~ 0.3, adopting a standard deviation of 07 ~ 0.04/arcsec for 
the flexion gives similar signal-to-noise ratios for the shear (1/3) 
and the flexion (1/2). Moreover, combining shear and flexion in- 
formation in the analysis can provide tighter constraints of the 
halo properties lEr et al.ll2012l) . 

A high number density of background galaxies is crucially 
important. However, if the two types of flexion can be mea- 
sured with sufficient accuracy, we shall already be able to place 
tight constraints on the halo ellipticity with current surveys, e.g. 
the HST Cosmic Evolution Survey (COSMOS) and the Canada- 
France-Hawaii Telescope Legacy Survey. We thus look forward 
to future projects with high precision images such as the Dark 
Energy Survey and the Euclid mission. 
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